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Abstract 

We present preliminary results on the Dalitz plot analysis of — > K~^7r~7:'^ decays. The data 
sample comprises 213 million T{4S) — > BB decays collected with the BABAR detector at the PEP- 
II asymmetric-energy B Factory at SLAC. We report measurements of the inclusive branching 
fraction, quasi-two-body fractions and CP-violating charge asymmetries for intermediate states 
including K*{892)+tt- and p{770)-K+. Observations of 5° decays to the Ktt S -wave intermediate 
states, i^:^(1430)+7r~ and i^o*(1430)°7r°, are reported. Evidence of the decay B° ir*(892)°7r° is 
seen. We set upper limits at 90% confidence level on branching fractions of the nonresonant and 
other less significant intermediate states. 
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1 INTRODUCTION 



Charmless three-body B decays significantly broaden the understanding of B meson decay mech- 
anisms and provide additional possibilities for direct CP violation searches. The B^ K^tt~tt^ 
decay is known to have contributions from the charmless intermediate B^ i<'*(892)"'"7r~ and 
S° p{770)~K+ decays. Although not yet observed, B^ Er*(892)°7r° could also contribute, 
as well as K^{U30), i^|(1430), K*(1680), /o(1450), and p(1700). The 5° ^ K+tt-tt^ decay can 
also occur via long-lived charmed intermediate states, i.e. , B^ — > D^tt^, and the doubly-Cabibbo- 
suppressed B^ — > D~ , D~ — > 7r~7r'^. Intermediate states with broad resonances are interfering 
and therefore a full amplitude analysis is required to extract the amplitudes and relative phases. 
Quasi- two-body branching fractions and CP-violating charge asymmetries for intermediate states 
are also measured. The results can be used to probe the weak phase 7 = arg [—yudyub/^cdYd^ 
the Unitarity Triangle ^ |2j . 

The charmless intermediate states are largely dominated by a gluonic penguin diagram (6 — > uus 
or dds), while the electroweak penguin contribution may also be sizeable jS]. The Dalitz plot analysis 
provides a sensitive probe on the penguin contributions to B decays to a charmless vector meson 
and a charmless pseudoscalar meson. 

Measurements of some of the branching fractions and charge asymmetries have been carried 
out by CLEO BABAR 5 and Belle [S] in the quasi-two-body approximation. In this report, 
we present a Dalitz plot analysis that takes into account interference between individual modes, 
using a Maximum Likelihood (ML) approach. Unless stated otherwise, charge conjugation is always 
implied throughout this paper. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy e'^e~ storage ring. The sample consists of 213 million BB pairs, corresponding to an 
integrated luminosity of 193.2 fb~^ collected at the T(4S') resonance ( "on-resonance" ) , and an 
integrated luminosity of 16 fb~^ collected about 40 MeV below the T(45') ("off-resonance"). 

The BABAR detector is described in detail elsewhere (Tj. A five-layer double-sided silicon vertex 
tracker (SVT) and a 40-layer drift chamber (DCH) are used to detect charged particles and measure 
their momentum as well as ionization energy loss [dE/dx). Charged hadrons are identified with 
a detector of internally reflected Cherenkov light (DIRC) and ionization in the tracking detectors. 
Photons, neutral hadrons, and electrons are detected in a CsI(Tl) calorimeter (EMC), while muons 
are identified in the magnetic flux return system (IFR). 

3 ANALYSIS METHOD 

3.1 Event Selection 

Each signal B candidate is reconstructed from , vr", and vr'^ candidates. Charged tracks must 
have at least 12 hits in the DCH and a transverse momentum larger than 100 MeV/c. Charged 
tracks identified as electrons, muons or protons are rejected. We also require that the charged kaon 
candidate must be identified as a kaon and that the charged pion candidate must not be identified 
as a kaon. The vr" candidate must have a mass that satisfies 0.11 < 772(77) < 0-16 GeV/ (?, and each 
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photon is required to have an energy greater than 50 MeV in the laboratory frame and to exhibit 
a lateral profile of energy deposition in the EMC consistent with an electromagnetic shower. 

Two kinematic variables, IS.E and mESi allow the discrimination of signal B decays from random 
combinations of tracks and vr*^ candidates. The energy difference, A.E, is the difference between 
the e"'"e~ center-of-mass (CM) energy of the B candidate and -v/s/2, where ^/s is the total CM 

energy. The beam-energy-substituted mass, mESj is defined by \J {s/2 + pj • /Ef — p^, where 
the B momentum, pe, and the four-momentum of the initial e'^e~ state {Ei, pj) are measured in 
the laboratory frame. For the signal B candidate, we require that mES > 5.27GeV/c2. The /^E 
resolution depends on the vr*^ energy and therefore varies across the Dalitz plot. To account for this 
effect, we introduce a transformed quantity: 

AP AP ' 

where /^E^ax = 0.08 - 0.0014 • m(K+7r"), AS^m = -0.14 + 0.0038 • m(K+7r"), ah in units of 
GeV. We require \l\E'\ < 1. 

Continuum e'*"e~ QQ {q = u,d,s,c) events are the dominant background. To enhance dis- 
crimination between signal and continuum, we utilize a neural network (NN) which combines four 
discriminating variables: the cosine of the angle between the B momentum and the e"*" beam direc- 
tion in the CM frame, the cosine of the angle between the B thrust axis and the beam direction in 
the CM frame, and two event-shape variables defined as sums over all particles i of pi x \ cos^jl*^, 
where n = or 2, and 9i is the angle between the momentum of particle i and the B thrust axis. 
The NN weights the discriminating variables differently, according to training on off-resonance data 
and simulated signal events. We cut on the NN output to suppress about 95% of the continuum 
background while retaining 62% of the signal. 

The fraction of events that have more than one candidate satisfying the selection depends on 
the kinematical distributions of the B decay products, and is therefore different for each mode. 
When more than one vr*^ candidate is present in the event, we choose the candidate with the 
reconstructed ir^ mass closest to the nominal value of 0.1349 GeV/c^. If more than one charged 
track combination satisfies the selection, we randomly choose one B candidate. An event is classified 
as misreconstructed signal if it contains a B which decays to the signal mode, but one or more 
reconstructed particles are not actually from the decay of that B. Misreconstructed signal is called 
self cross feed (SCF). This misreconstruction is primarily due to the presence of low momentum 
pions in the B decays and thus the SCF fraction varies across the Dalitz plot. 

There are 7220 events selected from the data sample. Tabled summarizes the signal efficiencies 
(including the misreconstructed events) and the misreconstruction fractions (with respect to the 
total signal) estimated from simulation with non-interfering amplitudes. 

We use simulated events to study the background from other B decays (.B-background) which 
includes both charmed (6 c) and charmless decays. In the selected B^ — > K~^7r~7r^ sample we 
expect 503 it 151 b ^ c background events. The decays B^ — > D tt^^D — > K^tt~ and B^ 
D~ , ir^TT^ have the same final state as the charmless intermediate states of our interest, 

but they do not interfer with the charmless intermediate states. Both are technically treated as 
S-backgrounds, but separately from other charmed background events. Branching fractions of 
S° ^ TTtt^ and 5° ^ D-K+ have been already measured, albeit with large uncertainties. In 
the present analysis, while we still do not have enough sensitivity to observe the doubly-Cabibbo 
suppressed B^ D^K^, ir^n^, we measure the branching fraction of the B^ D^vr^ 

decay. The selection efficiency for this mode is (16.7 ib 0.1)% with a misreconstruction fraction of 
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Table 1: Signal efficiencies (misreconstructed events (SCF) included), fractions of misreconstructed 
events, for the different intermediate signal modes. The errors given are statistical only. 



Decay mode 




/scf(%) 


i^*(892)+7r- 


14.9 ± 0.1 


23.6 ± 0.3 


p-K+ 


16.7 ± 0.0 


22.2 it 0.1 


i^*(892)07r° 


16.8 ± 0.1 


8.7 ± 0.2 


nonresonant 


16.3 ± 0.0 


3.8 ± 0.0 


iro*(1430)+7r- 


15.5 ± 0.1 


9.0 ± 0.2 


i^^(1430)°7r° 


17.4 ± 0.1 


4.2 ± 0.2 



(2.6 lb 0.1)%. We expect 230 ± 20 charmless 5-background events in the data sample. The major 
charmless B-background modes include 5° ii'*(892)°7, K*(1430)°7, K+Tr', B+ i^+7r°, 

p'^7r~, B^ P'^P~, and 5+ p'^ ■ Other contributing modes include 5+ ir*(892)+7r°, 
5+ f^K+, B+ ii:Q(1430)°7r+, and B+ -> r]'K+. The branching fractions of unmeasured 
decay channels are estimated within conservative error ranges. Charmless -B-background modes 
are grouped into six classes with similar kinematic and topological properties. Two additional 
classes account for the inclusive neutral and charged 6 — c decays. In total, ten classes for B 
backgrounds are defined. 

3.2 Dalitz Plot 

For the decay B^ — > K~^Tr~Tr^, with four momentum p^o, p+, p- and po respectively, the differential 
decay width with respect to the invariant-mass-squared variables (i.e., the Dalitz plot variables 
(DP)) reads 

irlB^ K^n-.') = ^t^-^-^ , (2) 

where Abo—^k+-k--!t° is the Lorentz-invariant amplitude of the three-body decay and 

s+- = {p++P-f , s-o = {p-+pof . (3) 

The amplitude Ab"^k+tt-tt" contains all the underlying dynamics of the B^ K~^tt~tt^ decay. 
In general, it is the coherent sum of one nonresonant (NR) term, ^nr, which is assumed to be 
constant in the Dalitz plane, and of several resonant amplitudes, having spin J, magnitude Oj, 
and phase 4>i- 

Abo^k+'k-kO {s+- , s-o) = ONRe^'^'^^^NR + ^ aic'^^ ■ A(s+- , s-o) . (4) 

i 

A similar expression is implied for B^ decays, with onr, ^nr; being assumed different 

from Onr, <;i!'nR) and (pi, respectively. 

The resonant amplitude Ai is written as a product of four terms 

A(DP) = '^FB,i ■ -^Fiis) ■ •^K(DP) • '^Fr^s) , (5) 

where "^Fb^i is an irrelevant constant form factor for the B° decay and is absorbed by normaliza- 
tion, '^Fi{s) = '^F{Rq{s)) /'^F{Rq{mf)) is the ratio of Blatt-Weisskopf damping factors (see below). 
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is a kinematic function (see below), and Fji^i{s) is a dynamic function describing the 
resonance (see below). 

For a resonance decays to particles a and b, with invariant mass ^/s, the momentum of the 
particles a and b in the resonance rest frame, q, is given by 

q{s) = ^ (l - ^^^I^^^±^Y (^1 - (jlh^^^Y , (6) 

where rua and m?, are the masses of particle a and b, respectively. 

The functions '^F{Rq{s)) are the nuclear Blatt-Weisskopf damping factors jH], given by 

Of = 1 ^F = ^ ^F = ^ (7) 

where R is the range parameter. The Blatt-Weisskopf damping factors are studied only for sys- 
tematic uncertainty evaluation. 

The spin-dependent function is taken to be equal to 1 for a spin-0 resonance, — 2|pd||pb| cos9 
for spin-1, and |(|pd||pb|)^(3cos^ — 1) for spin-2, where pd is the three-momentum of one of 
the resonance daughters and Pb is the three-momentum of the bachelor particle, both measured 
in the resonance rest frame, and 6 is the helicity angle of the resonance. For a resonance formed 
from K^ir^ (tt^tt^), the helicity angle is defined by the angle between the vr^ (vr^) in the resonance 
rest frame and the resonance flight direction in the rest frame. For a resonance formed from 
K~^TT~, the helicity angle is defined by the angle between the in the resonance rest frame and 
the resonance flight direction in the rest frame. 

Three parameterizations are considered to describe the decay dynamics, Fji^i{s). Parameters 
are taken from (2] unless stated otherwise. The relativistic Breit-Wigner parameterization with 
mass-dependent width is used for K*(892)+'0, K|(1430)+'°, and K*(1680)+'°. It is given by 

'Fn,{s) = 2 jr f ^ ' (8) 

s — + inii-^ L i[sj 

The s-dependent width is defined by 

= (^y'^' 'F\Rq{s)) 

'^s \q{mi) J JF^{Rq{m1)) ' 

where rrii is the mass of the resonance i, F-* = Fj(m?) its width. 

The Gounaris-Sakurai parameterization jTU] is used to parameterize p(770)^, /9(1450)^ and 
p(1700)^. For the iTvrS'-wave resonances, -fCg (1430)"''''', which are found to be dominant in the 
Kt: invariant mass range below 2 GeV/c^, an effective-range parameterization was suggested 
to describe the slowly increasing phase as a function of Kt: mass. We use the parameterization as 
used in the LASS experiment tl2j, being modified for B decays: 

^LA55/ ^ ^ , g2i<5B F?/ g(mf ) 

^' q{s) cot 6b - iq{s) mf - s - imiFiis) ' ^ ' 

where ^ ^ 

cot^B = — -- -h -rg(s), (11) 
aq[s) 2 
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nii = 1415±3MeV/c2, = 300±6MeV/c2, the scattering length a = 2.07±0.10 (GeV/c)-i, and 
the effective range r = 3.32 ± 0.34 ( GeV/c)~^. Our results for ifQ(1430) reported in this paper are 
not purely due to the resonant term, but to the K-k S-wave as a whole. 

Our nominal model includes the nonresonant contribution and five resonant intermediate states: 
p(770)-, i^*(892)+'0 and i^^(1430)+'°. Variations to this nominal model are used to estimate the 
model-dependent systematic uncertainty in the results. 

The quasi-two-body fractions and CP-violation charge asymmetries are extracted as follows: 

^ ^ / ( |a»e^'^' + iQje^^' ) ds+.ds-o ^^2) 

/ ( \Ab0^k+^-^0 |2 + |^gO^^_^+^o |2 ) ds+-ds-Q ' 

where Oj and aj are the fitted magnitudes, (pi and cp^ are the fitted relative phases, for the interme- 
diate state i. The extraction is also valid for the nonresonant. Due to interference, the sum of the 
fractions will in general not add up to unity. 

The inclusive branching fraction, B^"''^'", and the quasi-two-body branching fraction for an in- 
termediate state i, Bi, are given by: 

where Nf^ig is the total signal observed in the data, e is the signal efficiency averaged over the 
Dalitz plot, N^-^ is the total number of BB pairs produced, and fi is the fitted fraction for the 
intermediate state i, respectively. 



3.3 Maximum Likelihood Fit 

The selected on-resonance data sample is assumed to consist of signal, continuum-background and 
5-background components. We use the variables tties, Ai5^', and the Dalitz plot to discriminate 
signal from background. The signal probability density function (PDF) contains two parts cor- 
responding respectively to signal events that are correctly reconstructed (TM) and signal events 
that are misreconstructed (SCF). A PDF is introduced to describe the dominant continuum back- 
ground. For S-related backgrounds, we study the contributions with exclusive simulations and 
group similar contributions into the 10 B background classes, each of which has an individual PDF 
in the fit. 

The likelihood C for N events reads as 

, N ^ 

(1 ~ /sCF)^s«ff-TM,i + /sCF^s«ff-SCF,i 

i=l ^ 

+ (1 + 9K,i-4gg) Vqq,i 

classes -| 

3=1 



c = e-^ n 
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Nf 

classes 



where, 



E ^Ko\{^ + <lK,A^Bo)Ko,} > (15) 



• N is the sum of all the yields involved, i.e., Nsig + Ngg + J2j=i"''' + Ej=i"" ^^o! 

• Naig is the total number of signal events in the data sample to be determined in the fit; 

• /scF is the fraction of the misreconstructed signal events averaged over the Dalitz plot, which 
is determined from the fit; 

• Vsig-TM,i and Vsig-scF,i are the products of discriminating variable PDFs, for the correctly 
reconstructed and the misreconstructed signal events, respectively; 

• Ngg is the number of the continuum events, which is determined in the fit; 

• qK,i is the kaon charge of the event; we use = 1 for and qk = —1 for B^; 

• Aqq parameterizes possible charge asymmetry in the continuum events due to detection, 
reconstruction or selection; it is free in the fit; 

• 'Pqq,! is the PDF for the continuum events; 

• -^dasses ("^dasses) the number of the charged (neutral) -B-related background classes consid- 
ered, equal to 6 (4); 

• iV^+ i^Bo) is the number of events in the charged (neutral) S-related background class j; it 
is fixed to the MC estimate unless stated otherwise; 

• A'^^ ( A'^o ) describes the charge asymmetry in the charged (neutral) B background of class 
j; this eventual charge asymmetry could come from physics, or detection effects; it is fixed in 
the fit; 

• 7^^+ j is the i3^-background PDF for class j; 

• "P^o j is the B'^-background PDF for class j; two dedicated classes are used for B^ — s- 
;dV ^ K+TT- and B° D-K+ ,D- tt-tt^. 

Each PDF Vx is a product of the PDFs of the four discriminating variables, i.e. Vx,i = 
PxM^s) ■ PxA^E') ■ Px,i(s+-, s-o). 

The TTiEs distribution for the correctly reconstructed signal events is parameterized using a 
Crystal Ball function T^, the mean and width being simultaneously determined from the fit to data. 
A sum of two Gaussians is used to describe the IS.E' distribution, the means, widths and fractions 
being parameterized as linear functions of the Dalitz plot variable m?(K^TT~), the intercepts and 
slopes of which are obtained from simulation, and are varied for systematic uncertainty studies. 
For the misreconstructed signal, the ttj-es and AE' distributions are both obtained from simulation. 

For the continuum events, the ttt-es distribution is described by an ARGUS shape function |14| 
and its AE' distribution is modeled by a linear function. The shape parameter ^ of the ARGUS 
distribution and the slope of the linear function are determined in a fit to the on-resonance data but 
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selected with m-Es > 5.23GeV/c^. This dedicated fit uses only uies Ai?' as the discriminating 
variables, and considers the signal, the continuum background, and the -B-backgrounds in a similar 
way as the nominal fit. 

The background class T^ir^ , K~^tt~ shares the parameterizations of ttt-es and AE' 

for the signal events. Its yield is determined in the fit to measure the branching fraction. 

The shapes of the other B background mss - A£" PDFs are obtained from MC simulation and 
parameterized using two-dimensional empirical shape- fitting techniques |15j . 

The Dalitz plot PDFs for TM and SCF signal events are 



A/tm 

PsCF,i = SihcF / K . I ,+,ol ^ ^^^^ 

■/VsCF 

where A/tm and A/scf are the normalization constants to be dynamically determined in the fit, Si 
and /scF,i are the efficiency and the misreconstruction fraction, respectively, that vary across the 
Dalitz plot and are obtained with phase-space MC simulation. The magnitudes and the relative 
phases in Eq. are fit parameters and are allowed to vary in the fit, except that we fix to 
1.0 and both the phases to for the p component, since this analysis is only sensitive to relative 
magnitudes and phases. The masses and widths for intermediate resonances are fixed in the fit and 
are varied by their uncertainties to estimate the systematic uncertainties in our results. Because the 
considered intermediate resonances are in general very broad compared to the Dalitz plot resolution, 
we ignore the resolution effect for the correctly reconstructed signal, while we apply a resolution 
function (matrix) for the sizable misreconstructed signal. 

The Dalitz plot PDF for the continuum events is obtained from the on-resonance sideband 
events selected with 5.20 < twes < 5.25GeV/c^, plus the off-resonance events selected with uies > 
5.20 GeV/c^, and is corrected for the expected S-background feed-through in the sideband of 5.20 < 
mES < 5.25GeV/c2. 

For the background class B^ D^n^ — > K~^7r~, we distinguish the correctly reconstructed 
events and the misreconstructed events. For the correctly reconstructed events, a single Gaussian 
is used for the reconstructed mass distribution, the mean and width being free in the fit, 
and an empirical 5t/i-order polynomial is used to account for the reconstruction efficiency effects 
for the reconstructed helicity distribution. For the misreconstructed events, a smoothed two- 
dimensional histogram is used to parameterize their contributions in the Dalitz plot, which are fixed 
to the MC expectation. For the other nine classes of S-background, we model their contributions 
in the Dalitz plot with smoothed histograms obtained from simulation. 

Multiple solutions occur in the fit for this Dalitz analysis. To study the effect, we randomize the 
initial values of the amplitudes and phases that are let free in the fit and then the fit is redone many 
times. We do observe that there are two sets of solutions that correspond to two log-likelihood 
(log C) values. The two log C values differ by 0.4 units. Parameters from the two solutions usually 
give very close values except for the phases of K*{892)^ and A'o(1430)'^. Our strategy is to choose 
the fit that gives the largest likelihood to obtain our nominal results. We also report the other 
solution for the phases of K*{892f and K^{U30f. 
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4 SYSTEMATIC STUDIES 



The systematic errors in the branching fractions and CP- violating charge asymmetries, due to 
assumptions about the decay dynamics, are referred to as "model systematic uncertainty". They 
are obtained by varying the resonance parameters within their uncertainties, and by adding less 
significant resonances to the nominal model. 

The other systematic errors in the branching fractions are obtained by adding in quadrature 
those evaluated for the signal yield, the systematic uncertainties in efficiencies of tracking, particle 
identification, vr'' reconstruction and those corresponding to other selection criteria. The other 
systematic errors in the Acp measurements are introduced by the uncertainties in the treatment 
of the B background and by possible charge biases of the detector. 

The systematic uncertainties in the signal yield are primarily due to the modeling of the signal, 
of the i?-background, and of the continuum in the ML fit. These are estimated by varying the 
parameters that are fixed in the nominal fit, within their uncertainties. The variations in the signal 
yield are added in quadrature. The parameters for the signal AE' PDF are allowed to vary in the 
fit to the data to estimate the corresponding systematic uncertainty. The expected yields from the 
background modes are varied according to the uncertainties in the measured or estimated branching 
fractions. Biases observed in toy MC fits are added in quadrature and assigned as a systematic 
uncertainty of the fit procedure, referred to as "fitting procedure" in Table El 

The basis for evaluating the systematic uncertainties on the cuts that are applied in the selection 
process is to study the differences in ttj-eS) A£" and the NN output between on-resonance data and 
Monte Carlo simulation. The corrections and uncertainties in the signal efficiencies are summarized 
in Table The parameters for the signal rriES ^-nd AE' PDFs are obtained from the fit to the 
data and are used to estimate the systematic errors due to the cuts on tties and AE', respectively. 
The difference between the data and Monte Carlo distributions of the NN output is extracted from 
fully-reconstructed — > K~^7r~TT^TT~^ decays via intermediate or D*, and is used to estimate 
the systematic error due to the cut on the NN output. 

Since B background modes may exhibit direct CP violation, the corresponding parameters (the 
A-g's in Eg. I15j) are varied by ±0.5 to conservatively estimate the uncertainties. 

Table 121 summarizes the various sources contributing to the systematic errors in the branching 
fractions. The dominant systematic errors are due to the vr^ resonstruction, the AE' cut, and the 
cut on the NN output. Table |31 summarizes the possible sources contributing to the systematic 
errors in the charge asymmetries. 

5 RESULTS 

The preliminary results are obtained with a fit to the data sample using our nominal model. We 
obtain 1230 it 74 signal events and 454 it 24 B^ D 7r^,D K~^7r~ events from the data sample. 
The results obtained with the nominal model are shown in Tabled The total sum of the fractions 
obtained with the nominal fit is (102.6 it 8.0)%. Table El shows the results for the less significant 
modes, obtained by adding in turn one mode to the nominal fit. The results obtained from these 
altered fits, for the modes in the nominal model, are compatible. 

The inclusive signal efficiency is estimated to be (16.6 it 0.1)% with the observed Dalitz plot 
structure. We measure the inclusive branching fraction to be: 

B{B^ K+TT^n^) = (34.9 ± 2.l{stat.) ± 3.9{syst.)) x 10"^ . (18) 
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Table 2: Breakdown of systematic errors for the branching fraction measurements. The error from 
each source varies shghtly for each intermediate state and the most conservative estimate is quoted. 





Charmless yield (events) 


D^-K^ yield (events) 


Continuum mss PDF 


4.3 




0.0 


Continuum /S.E' PDF 


2.0 




0.0 


Signal A^' PDF 


50.0 




12.4 


B backgrounds 


2.7 




11.2 


Fitting procedure 


22.7 




0.0 


Sub-total (relative) 


55.2 ( 4.5% ) 




16.7 ( 3.7% ) 




Efficiency and scaling systematics (relative) 


Tracking efficiency correction 




1.6% 




PID for tracks 




2.0% 




tt" reconstruction 




6.4% 




/S.E cut efficiency 




4.6% 




niES cut efficiency 




0.4% 




NN cut efficiency 




6.0% 








1.1% 




Sub-total 




10.3% 




Total systematic error 


11.2% 




10.9% 



Table 3: Breakdown of systematic errors for the Acp measurements. 





i^*(892)+ 


p(770)- 


iv:o*(1430)+ 


i^o*( 1430)0 


iC* (892)0 


NR 


Detector bias 
B backgrounds 


0.02 


0.04 


0.01 

0.06 


0.07 


0.07 


0.20 


Total systematic error 


0.02 


0.04 


0.06 


0.07 


0.07 


0.20 
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We find 

B{B^ ;D%°) = (3.3 ± 0.2{stat.) ± OA{syst.)) x 10"^ . (19) 

The signal significance for B{B^ i^*(892)''7r''), including statistical and systematic errors, is 
4.2(7. We set an upper limit of B{B^ K^tt^-k^ nonresonant) < 4.7 x 10^^ at 90% confidence 
level, with systematic errors being taken into account. 

There exists another solution with a log C 0.4 units lower than the maximum log C. This 
solution gives very compatible results except for the relative phases of iCQ(1430)'' and i^*(892)''. 

The Dalitz plot for the selected K^tt~tt^ events is shown in Fig.^ The tties and AE' projections 
are shown in Fig. [2j Projection plots of the invariant mass pairs are shown in Fig. El 
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Figure 1: Dalitz plot for the selected K~^7r tt^ events. The visible narrow band comes from the 
5° ^ decays. 



6 SUMMARY 



Due to lack of knowledge of the final state interactions in B ^ K^tt~'7t^ decays, this Dalitz analysis 
assumes uniform phase space for the nonresonant decay amplitude and utilizes parameterizations 
obtained from non-i?-meson experiments to model the intermediate resonant amplitudes. Non- 
unique parameterizations of the decay amplitude are covered by the systematic uncertainties. 

A comparison with other measurements is summarized in Table 1^1 

The measured branching fraction of B^ D^n^ is in good agreement with the current world 
average, (2.7 ± 0.8) x 10"^ 9;. 

The measured branching fractions and charge asymmetries of the decays B^ — > p{770)~ and 
B^ i^*(892)"'"7r~ are in agreement with the previously reported measurements [HISIE]- 

We report observations of the intermediate 5-wave decays, B^ — > (1430)"'"7r~ and B^ — > 
Ko(1430)°7r°, with the LASS parameterization [12]. We do not see any significant direct CP 
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Figure 2: tues (left) and AE' (right) projections for the selected K^tt^tt^ events. The data are the 
points with error bars. The fit result is the (top) solid line. The (bottom) dashed line represents 
the continuum background and the (middle) dotted line is the B background added on top of the 
continuum background. Note that the B backgound includes the B^ D^ir^, i^^^vr" decays. 



violation in these two decay modes. 

We observe the first evidence of B^ — > i^*(892)'^7r'' decay with a branching fraction in agreement 
with the reported upper limit p. We do not see any direct CP violation in this decay mode. 

We set upper limits on the branching fractions of the nonresonant decay and the other six less 
significant intermediate states. 
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(statistical only). The total fit result is the solid (red) line. The continuum background is the light 
shaded (yellow) and the i?-background is the dark shaded (blue) added on top of the continuum 
background. 
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Table 4: Results obtained with a fit to the data using the nominal model. The first errors in the 
table are statistical while the second and the third (when present) are systematic errors as discussed 
in Sec. ^ The last column indicates the likelihood change when the mode is omitted in turn in 
the nominal model. / and Ac p have been defined in Eq. 1121 while cj) and (j) have been defined in 
Eq. El For the other solution, ^i^:* (1430)0 = (106 ± 44)°, 4>k*{892)o = (-44 ± 48)°, and all the other 
parameters give very compatible results. 



Mode 


/(%) 


(j) / (j) (degree) 


Acp 


Alog£ 


i^*(892)+7r" 


10.41IJ ± 0.8 


138 ± 35 / 
174 ± 42 


-0.25 ±0.17 ±0.02 ±0.02 


-42.5 




24.6^2:9 ± 0.6 


0.0 (fixed) / 
0.0 (fixed) 


0.13io:l7 ±0.04 ±0.13 


-48.3 


iC^(1430)+7r- 


32.2 ±3.8 ±9.4 


115 ±34 / 
149 ± 35 


-0.07 ±0.12 ±0.06 ±0.05 


-92.9 


i^o*(1430)°7rO 


22.5 ±4.0 ±7.2 


-12 ± 40 / 
8 ±42 


-0.34 ±0.15 ±0.07 ±0.08 


-71.4 




5.8ti:^ ±0.6 


-160 ±41 / 
-144 ±41 


-O-Oliai ±0-07 ±0.11 


-23.8 


Nonresonant 


7.ll|^ ±0.1 


55 ± 28 / 
79 ±28 


-0.12^[!|^± 0.20 ±0.07 


-7.3 



Table 5: Results obtained for the less significant intermediate states when they are added to the 
nominal model in turn. The last column indicates the likelihood change. 



Mode 


Fraction (%) 


Alog£ 


i^2*(1430)+7r- 




7.3 


/^2*(1430)°7r° 


1 6+i-'^ 


3.4 


i^*(1680)+7r- 




5.8 


ivT* (1680)°7r° 




1.6 


^(1450)-K+ 


5.0i:? 


3.6 


p(1700)-K+ 




2.1 
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Table 6: Comparison between measurements of branching fractions (xlO^^) and CP- violating 
charge asymmetries. Note that Belle measured K"^, which is defined to have l.lGeV/c^ < 
m{K'^Tr~'^) < 1.6GeV/c^ and has a helicity angular distribution indicating a scalar behavior. 
Our preliminary results in this table are corrected for the secondary decay branching fractions 
except for those indicated by a f . Previously published results from BABAR on p^K^ (listed within 
parentheses) were obtained with a subset of 81 fb^^ of the data used for this analysis. 





CLEO 


Belle 


BABML 


B{B^ K+TT-TT^ inclusive) 
Acp{B^ —>■ i^+vr^vr" inclusive) 


< 40 
- 


36.6 ±4.2 ±3.0 
0.07 ±0.11 ±0.01 


34.9 ±2.1 ±3.9 
- 


B{B^ K*{892)+TT-) 
A (' rO , i<^* fQno\+^—\ 

A.qp[Jd — > J\ [OvZ) TT ) 


16 ±6 ±2 

n 9(7+0.33+0.19 
U-^'-'_o.34_o.08 


14.8 ± 4.6 ± 1.5 ±2.4 


10.9 ±2.3 ±1.5 
— U.zo it U.l ( it U.UZ it U.UZ 


K*f892)°7r°) 

A-^ \ J—' ' A. \. \ kJ i-^ j 1 \ 1 

^cp(5° ^i^*(892)V) 


< 3.6 


< 3.5 


30±09±05 

-o.oilj]:^^ ±0.07 ±0.11 


B{B^ p{770)-K+) 
Acp(-D p(770) K^) 


16±8±3 


15.1 ± 3.4 ± 1.5 ±2.1 

n OO+0.22+0.06 
^•^'^-0.23-0.02 


8.6 ± 1.4 ± 1.0 
(7.3 ±1.3 ±1.3) 

n 1 O+0.14 1 r\ r\A 1 r\ 1 O 

O.lolg.iy ^ ^-^4 ± U.13 
(0.18 ±0.12 ±0.02) 


BiB'^ ^ Ko*(1430)07rO)t 
AcpiB^ K^iimfiT^) 




n ,+1.6+0.5 
"•-^-1.5-0.6 


7.9 ± 1.5 ±2.7 
-0.34 ±0.15 ±0.07 ±0.08 


B{B^ Ko*(1430)+7r-)t 
^Cp(5°^^o(1430)+^-) 




5.1 ± 1.5lg;f 


11.2 ± 1.5 ±3.5 
-0.07 ±0.12 ±0.06 ±0.05 


B{B^ K+TT-TT^ NR) 




< 9.4 


< 4.6 


B{B'^ i^2*(1430)+7r-)t 






< 2.2 


B{B^ ^ KKUSOfiry 






< 1.2 


B{B^ ii:*(1680)+7r-)t 






< 2.5 


B{B^ ^ if*(1680)07rO)t 






< 1.3 


B{B^ /9(1450)-ii:+)t 






< 3.2 


B{B^ p{nm)-K+)^ 






< 1.7 
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